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We report on numerical studies into the interplay of disorder and electron-electron interactions within the 
integer quantum Hall regime, where the presence of a strong magnetic field and two-dimensional confine- 
ment of the electronic system profoundly affects thermodynamic and transport properties. We emphasise 
the behaviour of the electronic compressibility, the local density of states, and the Kubo conductivity. Our 
treatment of the electron-electron interactions relies on the Hartree-Fock approximation so as to achieve 
system sizes comparable to experimental situations. Our results clearly exhibit manifestations of various 
interaction-mediated features, such as non-linear screening, local charging, and (;-factor enhancement, im- 
plying the inadequacy of independent-particle models for comparison with experimental results. 
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1 Introduction 

The integer quantum Hall effect (IQHE) — observed in two-dimensional electron systems (2DES) subject 
to a strong perpendicular magnetic field [1] — has been well studied using single-particle arguments [2- 
11]. However, experiments on mesoscopic MOSFET devices have questioned the validity of such a simple 
single-particle picture. Measurements of the Hall conductance as a function of magnetic field B and gate 
voltage [12] exhibited regular patterns along integer filling factors. It was argued that these patterns should 
be attributed to Coulomb blockade effects. Similar patterns have been found recently also in measurements 
of the electronic compressibility k as a function of B and electron density ric in the IQHE [13] as well as 
the fractional quantum Hall effect (FQHE) [14]. From these measurements it turns out that deep in the 
localised regime between two Landau levels, stripes of constant width with particularly small k can be 
identified. These stripes consist of a collection of small-K lines, identifiable with localised states, and 
their number is independent of B. This is inconsistent with a single-particle picture where one expects a 
fan-diagram of lines emanating from (0, 0) in the {B, no)-plane. The authors of [13, 14] argue that their 
results may be explained qualitatively by non-linear screening of the impurity charge density at the Landau 
level band edges. Clearly, such screening effects — described within a Thomas-Fermi approach [13] — 
are beyond the scope of a non-interacting theory. This immediately raises a question on the status of the 
aforementioned universality of the QH transition [12, 15, 16], which was obtained largely within a single- 
particle approach [9, 17-23]. Indeed a recently developed theory of the IQHE for short-ranged disorder 
and Coulomb interactions shows that universality is retained due to J- invariance [24-26]. 

In the present paper, we quantitatively investigate the effects of Coulomb interactions on the compress- 
ibility in the {B. ric) -plane within a Hartree-Fock (HE) approach. HE accounts for Thomas-Fermi screen- 
ing effects while at the same time leading to a critical exponent v whose value is consistent with the results 
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of the non-interacting approaches [27-29]. We find that the observed charging lines in the compressibility 
can be well reproduced and that the width of each group of lines is well estimated by a force balance 
argument. Thus we can quantitatively explain the lines and stripes in the compressibility as a function of 
(-B, He), fully in support of the quaUtative picture proposed in [13, 30-32]. However, our results for the 
conductivity patterns of Ref. [12] do not seem to reproduce the observed stripes. This might be due to our 
stiU too small system sizes [33] or many-particle physics beyond the HF approximation [34]. 



2 The Hartree-Fock approach in Landau basis 

We consider a 2DES in the (a;, y)-plane subject to a perpendicular magnetic field B = Bez. The system 
can be described by a Hamiltonian of the form 

H^u.s = h^ + Vc= + ^ + Vi(r) + Vcir, r'), (D 

where a = ±1 is a spin degree of freedom, Vi is a smooth random potential modelling the effect of the 
electron-impurity interaction, Vc represents the electron-electron interaction term and m*, g*, and jUb are 
the effective electron mass, 9-factor, and Bohr magneton, respectively [35]. 
The electron-electron interaction potential has the form 

= il^^^ =E^c(q)exp[iq. (r - r')] , (2) 

with 



q 



Vcici) = / ■ 7 |. . (3) 

The parameter 7 will allow us to continually adjust the interaction strength; 7 = 1 corresponds to the bare 
Coulomb interaction. Using 7 < 1 merely serves as a numerical handle to the interaction strength and 
might be easily eliminated by rescaling of the disorder strength. = L'^ /2'itII is the number of states 
per Landau level, also referred to as the number of magnetic flux quanta. Choosing the vector potential in 
Landau gauge, A = BxBy, the kinetic part of the Hamiltonian is diagonal in the Landau functions [36] 



Vn,fc(r) = {r\nk) = —;==== exp 
V2"n!V7rZc-t' 



iky- ^ ^ 



2/2 



Hn ( 1 , (4) 



where n labels the Landau level index, k = 2Trj/L with j = 0, . . . , — 1 labelUng the momentum, 

Hn (x) is the nth Hermite polynomial, and Ic = y^h/eB the magnetic length. We use Ni Gaussian-type 
"impurities", randomly distributed at Vg, with random strengths Wg <E [—W, W], and a fixed width d, such 
that Vi{r) = Efii (ws/nd^) exp[-(r - r,)Vd2j ^ ^^(q) ^^^^^^ . j.) ^ith 

^i(q) = 2^ exp ( • r, 1 , (5) 

5=1 

where q^.y = ^irj /L and j = —N^, —N^ — 1, ... , iV^. The areal density of impurities therefore is given 
by ni = Ni/L"^. The Hamiltonian is represented in matrix form using the periodic Landau states \nk) and 
we have 



H^,fc;„',fc' = {nk\H^^^^\n'k') 

1 ag* n 

2 ^ ~rn. 



1 (7Cj 771 

" + o + ~i 1 f^'^c5n,n'5k,k' + V„ fc/ -I- fe, , (6) 
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with the cyclotron energy hoJc = heB/m*. The disorder matrix elements are given by 'Vn,k;n',k' = 
X^q {^)Sn,k:n' .k' (q), where mixing of Landau levels is included. The explicit form of the plane wave 
matrix elements Sn,k;n',k'{<l) = {nk\ exp(zq • r)|n'fc') and the Fock matrix FJ^ ^, can be found in 
Ref. [29]. The total energy E^ot in terms of the above matrices is given as 

i?tot=Tr(hD + iFD) = i^ ^ {2KMn',k-+KMn',k-)^l',k';n,k (7) 

<7 n^k'^n' ^k' 

where the density matrix is D„ b = .fa^n.a^m*h ^^^^ /" Fermi function and ^ the expan- 
sion coefficients for the Landau basis [29] . The sums n and n' run over the number of Landau levels, A'll, 
taken into account. A variational minimisation of (4'|i?2DEs|^) with respect to the coefficients yields 
the Hartree-Fock-Roothaan equation [37]. The simplest scheme for solving this self-consistent eigenvalue 
problem is the Roothaan algorithm [37, 38]. However, convergence of the Roothaan algorithm is rather 
poor. In most cases it runs into an oscillating limit cycle [35]. This Umit cycle can be avoided be minimis- 
ing a penalised energy functional 

-E'(Dold, Dncw) = -E(Dold, Dnew) +fe||Dold — Dnew|P , (8) 

instead of the actual HF energy functional where "old" and "new" denote two consecutive states in the 
HF-selfconsistency cycle. This then leads to the Level -Shifting algorithm (LSA) [39]. However, while 
LSA avoids limit cycles, it can and does lead to unphysical solutions. A further improvement is based on 
the Optimal-Damping algorithm (ODA) [38]. The iteration is carried out just as in the Roothaan algorithm, 
only that the new density matrix is a mixture of the old and the new density matrix, i.e. 

D = (1 - A)Doid + AD 

new ; (9) 

with a damping parameter A which is chosen optimally according to the direction of steepest descent in 
the total HF energy [35]. As it turns out the performance also depends strongly on the interaction strength 
and the position of the Fermi level. For some filhng factors and choices of parameters, we might find fast 
convergence with one of the algorithms described above. However, over the whole range of filUng factors 
only a combination of ODA and LSA can guarantee convergence in any case. In Fig. 1 we have depicted 
the convergence behaviour for a typical HF run at four different filling factors. Note that we have plotted 
the convergence as defined by Eq. (8), rather than the total HF energy, which is also sometimes used as 
a measure. In comparison to the energy which acquires a finite value at convergence, the error estimate 
I |Dnew — Doidl I goes to zero and its convergence is therefore easier to interpret. 



3 Wave function and charge distributions 

3 . 1 The local density of states 

The local density of states (LDOS) is a very interesting quantity which is experimentally accessible via 
spatial scanning tunneling spectroscopy (STS) [40^2] . The differential tunneling current between the STS 
tip and the sample is thereby proportional to the density of existing states at a certain energy. By definition 
the LDOS is the tunneling density of states (TDOS) weighted with the charge density at a position r, 
defined by 

U^OS{E,v)=Y,Wa{v)\H{E-E^), (10) 

a 

such that / (frlDOSiE.v) = J2a ^i^ " ^a) = p{E), where p{E) is the TDOS. In practice, the 5- 
function is broadened by temperature as well as an AC modulation voltage applied when using a lock-in 
technique for noise reduction. The energy window can then be described by a semi-ellipse around E with 
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Fig. 1 Dependence of the convergence on the number of HF iterations for the three HP algorithms. The four panels 
correspond to 4 different filling factors. 




Fig. 2 Local density of states for a non-interacting Hall system as a function of (a;, y) position as well as energy E 
for a system of size L = 300nm. The filling factor « 3 and B = 5T. Large LDOS values are indicated by white. 

the broadening AS, which is usually of the order of ImeV [41 J. Thus, the LDOS can be regarded as the 
charge density in a small energy interval and allows to study the charge distribution in the 3 dimensions (2 
spatial and energy). 
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Fig. 3 Local density of states as a function of x position as well as energy E for a system of size L = 300nm at 
constant y and B = ST. The filling factor z/ = (left) and v = 0.5 (right). Large LDOS values are indicated by white. 

In Figs. 2 and 3 we show results for a model of an n-type InSb structure of three-dimensionally dis- 
tributed donor atoms with a 2DES sitting on top of the cleaved surface [42]. The first 2 spin-split Landau 
levels are shown for Fig. 3. At filling factor v = 0, electrons are tunnehng from the tip into unoccupied 
states of the sample and one thus probes non-interacting physics. A spin- splitting equal to the Zeeman 
energy can be seen which is of equal magnitude in each level. In the right panel of Fig. 3, we show the 
situation at v = 0.5, i.e. when the Fermi energy lies right in the centre of the lowest Landau level. We 
observe that the high-LDOS features of this level appear to be pushed away from the Fermi energy, result- 
ing in the onset of the famous Coulomb gap ("Efros-Shklovskii gap", see [43]) in the TDOS. We would 
hke to note that the gap in the TDOS shall not be confused with the gaps appearing in the quasi-particle 
spectrum. In fact, it has been shown that despite the Unear Coulomb gap in the TDOS, the thermodynamic 
density of states remains finite at the Fermi level [29,44]. Furthermore, the fluctuations in this level are 
much stronger than in the other three levels. We interpret this as a signature of the screening mechanism, 
i.e. the electrons near the Fermi level screen the disorder such that the states in the higher levels only see a 
weakened disorder potential. We also see that the spin-splitting in the first level is enhanced as compared 
to the non-interacting case. Similarly, the spin-sphtting between levels 3 and 4 is also larger. These effects 
seems to be related to the direct and indirect spin splitting enhancements for instance recently reported in 
Ref. [33]. 



3.2 The charge distribution 
The spatial distribution of the total electronic density 

M 

"W=EEl^2(r)|' = i-'E E T.^nMn',k'Sn,k;n',k'{ci)eM-mr) (11) 
cr a=l cr n,k,n',k' q 

is also readily calculated in our approach. It details the screening mechanism by providing direct insight 
into the interplay of disorder and interaction. Let us start at the QH transition. Fig. 4 depicts the critical 
charge density at = 1/2 for a non-interacting system in units of no. The contour lines show the impurity 
potential Vi(r) where the critical energy Vi{r) ~ ep is highlighted by a thick line. The charge density 
evidently behaves according to the semiclassical approximation [21] and follows the equipotential Unes of 
Vi (r) . For the interacting case, however, we expect Thomas-Fermi screening theory to apply [45^8] . This 
approximation is appropriate for an impurity potential smooth on the scale of the magnetic length as well 
as a sufficient separation of the Landau bands, characterised by the condition hwc/lc > \J (|Vli(r)P). 
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Fig. 4 Spatial distribution of non-interacting electron density n(r)/no at S = 4T, 7 = and v = 1/2 as indicated 
by the colour scale. Solid contour lines show the equipotential lines of the Vi(r). The thick solid lines corresponds to 
ep. 

The induced electrostatic potential of the charge density 

m = f d^r'^P^ (12) 

47reeo J \r' - r\ 

and the impurity potential Vi(r) form a screened potential V^cr(r) — Vi(r) + e<j>{r). Here, n accounts 
for the positive background. Since a flat screened potential is energetically most favourable, one expects 
to find 14cr(r) = ep for the case of perfect screening. However, since fluctuations of the charge density, 
n(r), are restricted between an empty and a full Landau level, i.e. (A^ll — l)"-o < n{r) < A^ll^-o, where 
no = is the density of a full Landau level, the screening is not always perfect but depends on 

flie fluctuations in the impurity potential as well as on the filling factor [45^7]. With the condition that 
Tiujc is large compared to the potential energy scale, we can restrict the discussion to a single Landau level, 
i.e. A'll = 1- The plane can be divided into fully electron or hole depleted, insulating regions — where 
n(r) = or n(r) = tiq, respectively — and metallic regions — where ri(r) lies in between. Depending on 
the filling factor, the extent of those regions varies. Close to the band edge, insulating regions dominate. 
Screening is highly non-linear and transport virtually impossible. On the other hand, if disorder is weak 
enough, fliere exists a finite range of filling factors in the centre of each band where metallic regions 
cover most of the sample, percolate and render the whole system metallic. The disorder is effectively 
screened and transport greatly enhanced. In that case, the charge density nscr(r) can be obtained by 
Fourier transforming the screened potential. In 3D, this simply leads to the Laplace equation. For 2D [49], 
however, one obtains riscr(q) ~ ^(2eeo/e^)|q|^i(q) + ^^"-o'^g.OiWhere the |q| = only couples to the 
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homogeneous, positive background and thus does not contribute to screening of the impurity potential. In 
other words, in our model only the fluctuations 6n{r) are essential for screening. Hence, in 2D, a perfectly 
screening charge density would obey 

47reeo f ^ ,A2DV[(r') 

"■scr r = — / d^r'— — + lyriQ. (13) 

J |r — r'l 

Clearly, the actual charge density is expected to deviate from nscr(r) for several reasons. Firstly, the 
fluctuations of n(r) are restricted as discussed above. Secondly, (13) is valid for the Hartree case only. 
Taking the Fock contribution into account will introduce short wavelength fluctuations due to the tendency 
for crystallization. Fig. 5 shows results for the charge density of interacting electrons aiv — 0.5. Broken 
lines indicate the regions where nscr(r) exceeds the range for (5ri(r) either below or above, i.e. areas that 
cannot be screened at all and thus exhibit insulating behaviour Otherwise, we find the charge density to 




Fig. 5 Spatial distribution of HF-interacting electron density n(r)/no at _B = 4T, 7 = 0.3 and v = 1/2 as indicated 
by the colour scale. Contour lines show equidensity lines of (13). The broken lines indicate unscreenable (insulating) 
regions. The thick solid line shows nscr(i') = Tiscr = 0. 

follow riser (r) very closely. In this regime, the density is well described by (13) and the screening is very 
effective. Metallic regions dominate over insulating ones and transport is expected to be good. 

We hence have established that interactions at the level of HF already give rise to substantial differences 
in the LDOS and local charge density distributions. Thus one might wonder why theoretical models based 
on non-interacting charges work at all. However, comparing Figs. 4 and 5 side by side, we see that the 
charge density, following the percolating equipotential at half filling for the non-interacting Fig. 4, retains 
many of the properties of the HF-interacting charge density line given by riser (r) = riser — 0. Thus 
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non-interacting percolation-type arguments for network models [9, 11, 23] should indeed provide a good 
starting point for a qualitative description of the QH transition. 

4 Compressibilities and Coulomb blockade 

In order to understand the experimental results [13, 14] we need to focus on the physics of the 2DES 
near the Landau band edges. As opposed to the band centre, in this regime the measured compressibiUty 
exhibits very strong fluctuations as a function of electron density which is due to sharp jumps in the 
chemical potential upon varying the carrier density. A clearer picture of what the microscopic situation 
is can be obtained by a spatially resolved scan of k~^. UnUke the fixed-position {B, ne)-scans, the SET 
tip will thereby pick up the spatial variations of the chemical potential change. In the measurements of 
Hani and coworkers [13], spatially localised jumps of the chemical potential are clearly identifiable which 
form bent segments within the {x, ne)-plane. These jumps are fingerprints of local charge accumulation 
and allow to conclude that quantum-dot-like structures of electrons (near an empty Landau band) or holes 
(near a full Landau band) must be formed within the sample. The multitude of lines in the compressibility 
originating from the dots are indications of Coulomb blockade effects and provide information on the 
position, charging as well as spatial extent of the dots. Since the charge is confined to a small region 
within the sample. Coulomb repulsion amongst the particles requires a comparably large amount of energy 
for adding or removing an electron to or from the dot. An increase in the number of electrons is thus 
accompanied by equidistant jumps in the chemical potential. The bending of the compressibility patterns 
is a feature of the SET tip bias which acts as an additional potential at the tip site. As the tip scans across 
a particular dot, it affects its charging condition as a function of distance to the centre which results in the 
observed arc-like distortion. 

Now we want to investigate whether our model is capable of producing dot spectra as found in the 
experiment. The change of the local chemical potential with respect to the electron density can be computed 
for our system by noting that the local chemical potential is the functional derivative of the total energy 
with respect to the spatial electron density 



n] S [ [ ri'^ ri'^ ,?T.(r)n(r ) f j2„/ '^('' ) 

(5n(r) 5n{r) 



^(r) = ""77^ = — ^ / / d'vd'v' ^ ' \' = / d^r'^^ . (14) 



Hence, we confirm the anticipated results that the chemical potential is simply the electrostatic potential 
due to the electrons in the 2DES. In terms of our previous definitions, the electrostatic potential of the 
2DES reads 

0(r) = L-2^ J2 ^Vb(q)D^_,.„,_,,5„.fc;„,,fc,(q)exp(-zqr) (15) 

<7 n,k,n' ,k' q 

and the local inverse electronic compressibiUty can be evaluated as 

K-\r) a /x^+^(r) - ^i^{r) oc .^^+^(r) - ,^^(r). (16) 

We have calculated K~^(r) for the lowest two Landau levels for a system of size L = SOOnm without spin 
as a function of position and carrier density at magnetic field B = 2T. The SET tip potential, VTip(r), 
has been approximated by a Gaussian function VTip(r) = WTip/T^d'^ipCxp[-~{r — RTip)^/d^ip] with 
■"^Tip = 6meV and dxip = 5nm. The sample is then "probed" along the a;-axis at Rxip = {i/50)L, where 
i = 1, ... ,50. For each position of the tip the total potential has to be evaluated and a complete HE run is 
carried out since the tip affects the electron density at each position differently. The results are presented 
in Fig. 6. We find remarkable agreement with the experimental results. Distinct charging patterns appear 
at the Landau band edges, whereas in the centre of the bands the features are much weaker. Thus, we 
have shown that our model is able to exhibit Coulomb blockade patterns. We consequently also expect the 
{B, ne)-calculations to exhibit charging patterns similar to the experimental results and indeed in Ref. [29] 
we show good agreement between our model and the results of Ref. [13, 14]. 
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Fig. 6 Inverse electronic compressibility for a HF-interacting system of size L = 500nm as a function of position 
X and electron density rie in the two lowest Landau levels without spin. The y coordinate has been fixed at 250nm. The 
patterns correspond to jumps in the chemical potential due to charging of local charge puddles (see text for details). 



5 Interaction effects in the Hall conductivity 

The compressibility patterns found in the experiments discussed above were preceded by transport ex- 
periments that revealed patters in the conductance, also clearly interaction mediated. Cobden et al. [12] 
reported on Hall conductance measurements on mesoscopic MOSFET devices in the (Vg, i3)-plane, where 
Vg is the gate voltage, that displayed a contradictive behaviour to the expectations deduced from a single- 
particle model. The gate voltage Vg can be Unked to the electron number density as drie/dVg = C/e, 
where C is a material specific constant. A Une in the {Vg, -B)-plane originating from V^(ne = 0) can thus 
be labeled by a certain filling factor v. 

The fluctuations around the plateau transitions found in Ref. [12] at half integer filling factor are appar- 
ently correlated over a large interval of the magnetic field. Surprisingly, however, the lines formed by the 
conductance extrema in the {Vg, i?)-plane evidently ahgn with integer filling factors. This behaviour con- 
flicts with what is expected from a single particle model [50], where extrema in the transport coefficients 
can be linked to resonances in the disorder potential. 

Cobden et al. [12] suggested a scenario based on dividing the charge density in the sample into com- 
pressible and incompressible regions. Wherever the density is free to fluctuate, i.e. where the local filling 
27r^^ne(r) lies well away from an integer value, the 2DES exhibits metallic behaviour. Wherever the den- 
sity fills a whole Land level, i.e. 27rZgne(r) is close to an integer, the density is incompressible and thus 
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insulating. In the spatial density profile, the metallic regions form puddles that are enclosed by insulat- 
ing boundaries where the density corresponds to a full Landau level. With electron-electron interactions 
present, transport through the sample is now strongly influenced by those metallic puddles that are always 
surrounded by an incompressible density strip. The conductance peaks can then be associated with the 
charging condition of the puddles and therefore with the shape of the puddles. It is reasonable to assume 
that a particular density profile does not change along lines of constant filling factor in the {Vg, i?)-plane, 
only the average density as rio ~ veB/h. Thus it is clear that along lines where v is an integer, the shape of 
the incompressible puddles (contours where 27r/^nc(r) is an integer) remains roughly constant, therefore 
its charging condition and also the conductance extrema. 

Using the Kubo formula for the conductivity [51,52], we have calculated the Hall conductance in the 
{B, ne)-plane for different disorder and interaction strengths. In Fig. 7 we show the lowest two spin 
levels of a system of HF- interacting electrons of size L = 300nm. The disorder strength is chosen as 
W/ (P = 2.5meV with an impurity range of d — 40nm. The spin levels are well separated by the exchange 
enhanced spin splitting. We observe stable plateaus around integer filling factors, which are indicated 
by the black broken lines. We have also indicated the boundary between the non-linear and the linear 




m 

Fig. 7 Hall conductivity, axy, for the lowest two spin levels of a HF-interacting system of size L = 300nm with 
disorder strength W/d^ = 2.5meV in the (_B, tic) -plane. The two spin levels are well separated by virtue of the 
exchange enhancement of the spin splitting. Integer filling factors are indicated by the black broken lines, whereas the 
grey broken lines indicate boundaries between linear and non-linear screening [29] . 

screening regime by grey broken lines [29]. 

From our numerical results, we observe wide, stable plateaus at the integer filling factors. In contrast to 
a single-particle calculation [52] where the widths of the plateaus depend linearly on the magnetic field B, 
the widths of the plateaus remain constant as a function of B when electron-electron interactions are taken 



© 2003 WILEY- VCH Veikg GmbH & Co. KGaA, Wcinhcim 



pss header will be provided by the publisher 



13 



into account. Similarly to the compressibility calculations, the estimation of the cross-over between the 
non-linear and linear screening appears to describe the competition between disorder and interactions well. 
The constant width of the plateau is in agreement with the experimental finding of Ref. [12]. The expected 
aligrmient of conductance peaks along integer filling factor is, however, absent in our calculations. Instead, 
we observe seemingly random conductance jumps in the centre of the bands. We attribute this behaviour 
to the strong exchange correlation in this regime, which was also predominant in Ref. [29], where the 
compressibility became strongly negative in the centre of the bands. Therefore we conclude that exchange 
induced effects dominate over charging and Coulomb blockade effects. 

In Fig. 8 we show cross-sectional cuts in the {B, ne)-plane at -B = 4T and 6T. Results for 2 Landau 
levels of non-interacting (a,b), strongly disordered (c,d,e,f) and weakly disordered (g,h,i,k) systems are de- 
picted, corresponding to 2 spin-less Landau levels (c,d,g,h) or 1 spin-split Landau level (e,f,i,k). We have 
indicated the linear screening by the vertical grey shading. The plateau regions of localised electrons align 
well with the estimation formula of Ref. [29] and, to the accuracy of this simulation, their width is indeed 
independent of the magnetic field. The delocalised regime of the plateau-to-plateau transitions on the other 
hand increases strongly. This is in contrast to the single particle result, where both, the plateau and the tran- 
sition region increases with the magnetic field. Whereas the compressibiUty stripes are fingerprints of the 
non-linear screening, i.e. the insulating regime, the increasing width of the plateau transition demonstrates 
an interaction promoted delocalisation. We regard this behaviour as another important manifestation of the 
competing interplay between disorder and electron-electron interactions even in the integer QH effect. 

6 Conclusions 

We have investigated numerically how the interplay of electron-electron interactions and disorder affects 
transport and thermodynamic properties of electrons in the IQHE. We diagonalised the Hamiltonian for 
electrons confined to two dimensions and subject to a perpendicular magnetic field in the suitable basis 
of Landau functions and treated interactions in an effective, self-consistent HF mean-field approximation. 
Our results are qualitatively in very good agreement with recent experiments. More quantitative state- 
ments especially regarding spin- splitting enhancement might be limited by the unscreened HF approach 
used in this paper. Electron correlations beyond the exchange interaction may lead to a screening of the 
exchange term. However, a simple replacement of the exchange term by a screened one will destroy the 
self-interaction cancelation occuring among Hartee and Fock term, thereby introducing a new source of 
error. We believe that for the purpose of this investigation our approach describes the important physical 
effects well [44,53], especially for thermodynamic quantities where the system is given enough time to 
relax into equilibrium. The possibility to go beyond a simple electrostatic theory [54] by incorporating 
the exchange correlation effects, such as crystalisation and enhanced spin band separation, while at the 
same time allowing for the treatment of system sizes comparable to experiments, render the HF method a 
suitable choice. In that respect more complicated approaches such as density functional schemes [55] may 
not outweigh the advantages of the unscreened HF approach. 

We close with some final remarks. Whereas disorder is usually associated with a reduction of sig- 
nal quality, in the IQHE the concurrence of a magnetic field, reduced dimensionality, and disorder leads 
to a remarkable resihence of quantised transport. While the critical behaviour at the transition appears 
unaffected, localisation in the band tails is changed by interactions, yielding a significant change in the 
widths of the plateau regions. Although single-particle models can well describe many aspects of quantum 
Hall physics, only by taking interactions into account can the whole spectrum of experimentally observed 
features be satisfactorily understood. 
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Fig. 8 Hall conductivity, axy, as cross-sectional cuts for 5 = 4 (left column) and 6 (right column) and (a,b) two 
spin-less Landau levels with no Coulomb interaction, (c,d) and (g,h) a spin-split Landau level at W/ cf = 2.5meV and 
W/cP = 1.25meV, respectively. Panels (e,h) and (i,k) refer to two orbital, spin-less Landau levels at W/ = 2.5meV 
and W/d^ = 1 . 25mcV, respectively. The grey areas indicate the linear screening regime [29] which widens for weaker 

disorder (bottom two rows). ^ 2003 WILEY-VCHVerlag GmbH &Co. KGaA, Wetahem. 



